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1. Introduction 


While most projectiles in use today behave as rigid bodies while in flight, a notable number of 
projectiles are purposefully designed to carry a liquid payload. For example, smoke screens 
delivered by artillery rounds consist of a typical spin-stabilized shell containing a canister filled 
with white phosphorous { 1 - 4 ). Another example is new less-than-lethal projectiles having a 
concentric cylindrical cavity filled with liquid that delivers this payload to a target upon impact 
(5). Finally, some projectiles are designed to be general payload delivery shells, including 
delivery of medical supplies, such as intravenous fluid bags ( 6 ). 

There can be a significant difference in flight behavior between liquid- and solid-filled 
projectiles. The difference is caused by motion of the liquid inside the spinning projectile. This 
motion causes forces to act on the projectile body, which can prematurely terminate the flight by 
instability. Characteristics of this instability are sharp increases in angle of attack accompanied 
by large changes in spin rate { 7 - 10 ). At gun launch, the motion of the projectile causes the fluid 
to spin up in a time-dependent manner, but it subsequently achieves steady flow. 

Prediction of instability induced by a liquid payload installed in a coning projectile has been 
analyzed by computational fluid dynamic (CFD) simulations and with analytic linear fluid 
theories based on spatial-eigenvalue methods. The CFD solutions are generally best suited to 
liquids with a low Reynolds number (Re), while spatial-eigenvalue methods can encompass a 
much broader range of Re. A good survey of liquid-filled projectiles with a focus on flight 
instabilities is given in Sedney { 11 ). Generally, analytic theories are a composition of inviscid 
flow and viscous boundary-layer corrections { 12 - 14 ). The prediction of flight stability of a 
liquid-filled projectile has also been studied using tricyclic linear-projectile theory { 15 ). This 
analysis assumes the effect of a liquid payload is similar to the Magnus effect. Spectral analysis 
used to numerically compute liquid-fill induced moments has been incorporated into six-degree- 
of-freedom (6-DOF) simulations { 16 ). The angular motion of a low-Re liquid-filled projectile 
has also been simulated using a precomputed table of liquid moments, obtained from CFD, in 
6-DOF calculations { 17 ). 

The present report meshes a well-developed spatial-eigenvalue theory directly into a standard 
6-DOF projectile flight dynamic model. Hence, a well-established rotating liquid model 
calculates liquid moments at each time step of the 6-DOF trajectory numerical simulation 
without the need of interpolation to gain the effects of liquid payloads. This report further 
bridges the gap between the body of literature on effects of viscous liquid payloads on projectiles 
and 6-DOF projectile flight dynamic modeling. The report begins with a review of rigid 
projectile flight dynamic modeling along with a description of modeling a rotating liquid in a 
cylindrical cavity. The two models are subsequently integrated such that a projectile flight 
dynamic model with a liquid payload results. To highlight the utility of the methodology, the 
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flight dynamic model is exercised on an example shell. Comparisons are made between a liquid- 
filled projectile and a similar solid-filled projectile highlighting the predictive capability of the 
new model. 


2. Projectile Flight Dynamic Model With a Liquid Payload 


A typical 6-DOF rigid projectile model is employed to predict the dynamics of a projectile in 
flight. These equations assume a flat Earth. The well-known 6-DOF states comprise the three 
translational components describing the position of the projectile’s center of mass and the three 
Euler angles describing the orientation of the projectile with respect to the Earth. Figures 1 and 2 
provide a visualization of the degrees of freedom. 



Figure 1. Projectile position coordinate definitions. 
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Figure 2. Projectile orientation definitions. 


The equations of motion (18) derived in the no-roll frame are shown in equations 1-4. 
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Note that the following shorthand notation for trigonometric functions is used in equations 1—4: 
= sin(a), c a = cos(a), and t a = tan(a). 
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The force acting on the projectile in equation 3 comprises the weight force (W) and the 
aerodynamic force. The aerodynamic force is split into a standard force (A) and a Magnus force 
(M). The combination of forces is expressed in equation 5. 
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Equation 6 provides the expression for the weight force in the no-roll coordinate system. 
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Equation 7 provides the expression for the aerodynamic force in the no-roll coordinate system. 
This force acts upon the projectile at the aerodynamic center of pressure. 
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Equation 8 provides the expression for the Magnus force in the no-roll coordinate system. The 
Magnus force acts upon the projectile at the Magnus force center of pressure. 
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Equations 7 and 8 are based on Mach number-dependent coefficients and the total aerodynamic 
velocity given in equation 9. 

V = VfE + v“ + w" . (9) 

The moment acting on the projectile in equation 4 comprises the moment due to the standard 
aerodynamic force (A), the moment due to the Magnus aerodynamic force (M), the unsteady 
aerodynamic moment (EIA), and the liquid payload moment (L) as shown in equation 10. 
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The moment due to the aerodynamic force is expressed in equation 11. 
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The moment due to the Magnus force is expressed in equation 12. 
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The unsteady aerodynamic moments acting on the projectile are expressed in equation 13. 
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The coefficients used in this model are specific functions of the Mach number of the projectile. 
For fin-stabilized projectiles, Magnus force and moment are typically ignored since its effect is 
rather small for slowly rolling projectiles. 

The angular motion of a projectile is altered by the inertial waves propagating through the liquid 
payload. In turn, these waves act on the walls of the liquid cavity and generate a liquid moment. 
The impact of this moment can have a devastating impact on the projectile’s angular motion 
(19-21). For this report, the liquid moment is calculated from steady linear-projectile theory, 
which incorporates solutions of the linearized Navier-Stokes equations (18). In general, this 
liquid moment model is a linear combination of fast and slow mode liquid moment contributions. 
In section 2 of his report, Murphy (19) states that for a steady-state coning projectile with a 
liquid payload, “. . . only the fast mode motion is adversely affected by the liquid side moment.” 
For this reason, only the fast coning mode is used to calculate the liquid moment in this paper. 
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Figure 3 shows the orientation of body unit vectors (l B , J B , K B ) to velocity unit vectors 
(l v , J v , K v ). Appling the liquid theory (18) to figure 3 allows the liquid transverse and roll 
liquid moments to be written as shown in equation 14. 



Figure 3. Liquid transverse velocity orientation 
definitions. 


1 

tr 1 * 

r 

_i 


"l 

0 


0 

m l 

2 ~ 2 t-' 

= m L a p T 

0 

! 

V?* 

_n l ' 


0 




T 8 1 

f 

4 1 

( 

2 \ 

r 1 

^LRM 

"^LSM + ~ 


1 





2 



3 1 


) 


/v 2 +w 2 

V 





K=- 

-, P=- 







V 

w 






k 2 c 


LRM 


KC 

KC 


LSM 


LIM 


(14) 


The moment coefficients Clsm and Clim depend on the fast mode coning frequency T, fast mode 
undamping s, the cavity aspect ratio c/a, liquid Re, and the magnitude of axial spin (19). Values 
of the fast mode pair (T, s) are obtained from equation 15. 
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The liquid moment coefficients Clsm and Clim are calculated using the techniques given by 
Murphy (19). 


3. Results 


The example simulations given here consider liquid payloads in a typical projectile with launch 
properties given in table 1. 


Table 1. Launch properties of projectile. 


Physical Parameters 

Magnitude 

lx 

0.1157 slug ft 2 

Iy, Iz 

1.1972 slug ft 2 

m 

103.0 lb/ft 

D 

0.510 ft 

Launch velocity, V 

2460.0 ft s _1 

Launch axial spin rate, p 

1500.0 s- 1 


The liquid density p L is taken to be 3.5 times the density of liquid water, and the viscosity n is 
selected such that launch Re > 8 x 10 6 . This is sufficient to ensure that the Re remains large so 
boundary layer analysis (18) adequately governs the liquid physics for the entire trajectory. The 
liquid cavity is a cylinder with aspect ratio c/a = 3.75 completely filled with this low-viscosity 
liquid. Figure 4 shows the range of nondimensional coning frequencies, T, for a typical 
trajectory of this projectile housing a frozen (solid) liquid payload. 
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The Re for a given aspect ratio can cause large variations in the liquid moment when subjected to 
a range of nondimensional coning frequencies. Applying steady-state linear liquid theory to a 
payload configuration undergoing coning motion reveals in figure 5 the side moment coefficient 
C lsm with c/a = 3.75 for two numbers Re = 8 x 1 ()' : and 8 x 10 s . These results depict important 
C LSM behavior, where the peaks indicate a potential problem due to large liquid moments when 
the nondimensional coning rate T « 0.088 . Note the amplitude of the peak has a strong 
dependence on Re. Peak values of Clsm are nearly 2 orders of magnitude larger than the 
nonpeak values of Clsm- 



Figure 5. C L sm vs. coning frequency and Re, resonance T ~ 0.088. 

Figure 5 suggests that for a sufficiently large Re, a free-flight missile having nondimensional 
coning frequencies in the neighborhood of T = 0.088 will experience a significantly larger liquid 
moment than it does for frequencies outside this neighborhood. Such results are characteristic of 
the root cause of projectile instabilities due to liquid payloads. The aspect ratio, coupled with a 
large enough Re, forces inertial waves in the coning fluid to generate large C L sm coefficients 
(19). In general, increasing Re causes the liquid side moment to increase near T = 0.088, and 
this pronounced peak signifies a possible resonant frequency. This is one indicator of a liquid 
payload possibly causing the fast coning frequencies to change such that an instability occurs 
near T = 0.088. Note the range of frequencies given in figure 5 is close to the fast mode T 
frequencies for the frozen liquid found in figure 4. 

Next a series of results showing the effect a liquid payload has on projectile angular motion 
during flight is shown. Results presented here are for a launch Re = 8 x 10 6 , which is large 
enough so linear liquid theory with boundary layer corrections is valid throughout the entire 
trajectory. Figures 6 and 7 compare roll moment and spin rate results for a frozen liquid and a 
flowing liquid. These variables are selected because the rotational physics of a projectile causing 
a rapid decrease in spin rate is a strong indicator of flight instability caused by liquid payloads 
(18). 
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Figure 6. Projectile roll moment for frozen and flowing liquid payloads. 


The liquid payloads for the launch Re chosen in these examples has small effects on liquid roll 
moment, as shown in figure 6. As expected, figure 7 shows the projectile roll rates also change 
slightly compared to the frozen liquid payload. 



Figure 7. Projectile spin rate for frozen and flowing liquid payloads. 
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Evidently the increased magnitude of the projectile roll moment due to the flowing liquid is not 
sufficient to cause flight instability since the roll performance of the projectile changes very little. 

The liquid side moment coefficient, Clsm, is an important parameter when evaluating the impact 
a liquid payload has on projectile angular motion. Figure 8 shows the fast mode liquid moment 
coefficient as a function of T. This chart shows that Clsm changes rapidly for the region 
0.086 < T < 0.090 , and the magnitudes of Clsm exhibits the approximation C LSM ~ VRe (18) 
when the coning rate is in the neighborhood of T = 0.088. 



Figure 8. Liquid moment coefficients vs. nondimensional fast mode frequency. 

Figure 9 compares the angle of attack (AOA) for identical projectiles with payloads supporting 
inertial waves and frozen liquid payloads. 



Figure 9. Comparison of AOA of frozen and liquid payloads. 
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Apparently the inertial wave motion slightly increases the AOA for a liquid with launch Re of 
8.0 x 10 6 or 8.0 x 10 7 compared to the frozen liquid, but this increase is not enough to cause 
concern of flight instability. In fact, the results presented so far show no indication of flight 
instability even when the projectile has a coning rate near T = 0.088 for the potential problem in 
Clsm found in figures 5 and 8. 

In order to understand the dynamics of a liquid-filled projectile exhibiting angular instability, we 
consider figure 4 and exploit the peak value of Clsm near T = 0.088. Thus, to substantially 
increase the size of the liquid moment, we assume the payload comprises a hypothetical liquid 
with a decreased viscosity so that the launch Re is large enough to cause substantial increases in 
liquid side moment during flight. In particular, we want this liquid to generate large enough 
values of Clsm for T ~ 0.088 to alter the angular motion of the projectile. A sufficiently large 
decrease in projectile spin rate along with a large increase in AOA suggests the hypothetical 
liquid may cause premature termination of flight. Figure 10 shows the fast nondimensional T 
dependence of the side moment coefficient, Clsm, rapidly increasing near T = 0.088 for launch 
Re = 4 x 10 . This indicates a possible increase in the magnitude of projectile angular motion, 
which is consistent with the resonance configuration presented in figure 5. 



Figure 10. Comparison of projectile roll moment vs. payload launch Re. 


Figure 11 compares the AOA of the projectile with liquid payloads having three different launch 
Re’s. 
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The projectile with launch Re = 4 x 10 causes the fast frequency at t ~ 4.5 s, T ~ 0.088, to 
generate an AO A >50° for which the judgment was made to terminate the numerical integration. 
Figure 12 shows the corresponding projectile roll moment L rapidly increasing in magnitude due 
to the larger Re. Note that in the neighborhood of T = 0.088, the fast mode liquid motion is 
responsible for increasing the magnitude of L. 


Projectile Roll Moment 



Figure 12. Comparison of projectile roll moment vs. launch Re. 
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6 8 

Figure 13 shows the corresponding projectile roll rates for the launch Re = 8 x 10 and 4x10, 
which correlates with the size of liquid moment coefficient C L rm as shown in equation 15. The 

I I o 

increase in L for Re = 4 x 10 causes the rapid decrease in roll rate. 



o 

Figure 14 compares corresponding pitch rates of the frozen liquid and launch Re = 4 x 10‘ 
configurations. The rapid increase in q values near t ~ 4.5 s is caused by the fast mode liquid 
side moment. 
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A contrast in projectile yaw rates for frozen and liquid payload is given in figure 15. Again the 
rapid increase in yaw rate takes place near t ~ 4.5 s in the neighborhood of the fast mode 
resonance. 



These calculations show a significant difference in projectile angular motion when the liquid 

8 

payload has a launch Re = 4 x 10 compared to the same projectile with a frozen liquid. 


4. Conclusion 


An integrated trajectory simulation for a projectile with an internal liquid payload has been 
created and exercised. The liquid moments are calculated using linear liquid theory applied to a 
low-viscosity liquid undergoing two-mode steady state coning motion. High Re boundary layers 
approximate the liquid shear moment, yielding an average quasi-static liquid moment that is 
applied to a nonlinear 6-DOF time-dependent trajectory model. These quasi-static averages 
yield approximate predictions describing the motion of projectiles with low-viscosity liquid 
payloads filling a cylindrical cavity. The nonlinear 6-DOF motion of a projectile is often well 
approximated as the sum of fast and slow coning motions. Thus, under these conditions the 
linear liquid analysis can predict when a projectile exhibits flight instability due to a liquid 
payload by tracking the coning frequencies during the 6-DOF numerical integration process. 
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Calculations presented in this report show how both stable and unstable liquid payload 
configurations influence trajectory flight parameters. Unstable liquid payload configurations 
result from an improper combination of liquid Re’s, payload geometry, and quasi-static projectile 
coning mode rates. In such cases, projectile roll rate and axial moment decrease rapidly due to 
the quasi-static liquid moment. 
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List of Symbols, Abbreviations, and Acronyms 


6-DOF 

a 

AOA 

c 

C dd 

CFD 

I 


m 



Re 


t 

T 


six degrees of freedom 

radius of cylinder containing fluid (ft) 

angle of attack 

half length of the cylinder containing N subcylinders (ft) 

aero-coefficients 

computational fluid dynamics 

I.X 0 0 

projectile moment of inertia I = 1 0 I Y 0 [ (slug fr) 

° 0 I Y 

projectile mass (slug) 

liquid mass m L = 2 n a 2 c p (slug) 

projectile moment (ft lbf) 

projectile angular velocity (s ’ 1 ) 

displacement of d relative to mass center coordinates (ft) 

Reynolds number, a 2 p/// 
gyroscopic stability factor 
time (s) 

nondimensionalized coning frequency 
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projectile force (lbf) 

nondimensionalized growth rate per cycle 
projectile body Euler angles 
kinematic viscosity of liquid (ft s ) 

_o 

liquid mass density (slug ft ) 
variable in nonrolling system 
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